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D , In this paper we outline a new parton shower algorithm based on the Catani- 

Seymour dipole factorization. Our motivation is to have an algorithm which can 
naturally cooperate with the NLO calculations. 



£2 . 1. Introduction 



One often uses perturbation theory to produce predictions for the results 
of particle physics experiments in which the strong interaction is involved. 
In order to get useable predictions one has to calculate at least at next-to- 
leading order to avoid large uncertainties those come from the unphysical 
scale dependences. 

Unfortunately, standard NLO programs have significant flaws. One flaw 
is that the final states consist just of a few partons, while in nature final 
states consist of many hadrons. A worse flaw is that the weights are often 
very large positive numbers or very large negative numbers. 

There is another class of calculational tools, the shower Monte Carlo 
event generators, such as HERWIG [1] and PYTHIA [2]. These have the 
significant advantage that the objects in the final state consist of hadrons. 



* Contribution to the Proceedings of the Ringberg workshop on "New Trends in HERA 
Physics 2005" . This talk was given by Zoltan Nagy. 
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Furthermore, the weights are never large numbers. Finally, the programs 
have a lot of important structure of QCD built into them. For at least 
some cases like this, the shower Monte Carlo programs can provide a good 
approximation for the cross sections. The chief disadvantage of typical 
shower Monte Carlo event generators is that they are based on leading 
order perturbation theory for the basic hard process. This means that 
the predictions of the Monte Carlo programs usually have large theoretical 
uncertainty. From this reason we cannot use the Monte Carlo prediction 
as theoretical result in the experimental analysis. This rises an obvious 
question. Is it possible to build a Monte Carlo program that is predictable? 

The parton distribution function or the fragmentation function is a con- 
volution of a perturbative and a non-perturbative part. The perturbative 
part is the evolution kernel which is solution of the DGLAP equation. It 
is based on the collinear factorization property of QCD. Fixing the accu- 
racy of the calculation and the factorization scheme the evolution kernel is 
well defined. The conventional factorization scheme is the MS. The non- 
perturbative part cannot be calculated, fortunately it is universal and we 
can obtain it from fit to the data and use them for other processes. Of course 
if we change the factorization scheme we have to refit the non-perturbative 
part. 

What is the parton shower and the hadronization? We can think about 
that is a multi-hadron fragmentation function and the parton shower is its 
evolution kernel. The hadronization model is the non-perturbative part 
that we fit to the data. The perturbative part is based on the factorization 
theorem of the QCD and it is well defined up to the finite pieces. In other 
words, a different shower algorithm represent different shower scheme but 
upto the leading (LL) and next-to-leading (NLL) logarithmic accuracy they 
must be equivalent. If we have some free parameters other than the trivial 
scales in the parton shower algorithm then the different values of these 
parameters correspond to different shower schemes and the tunning of the 
hadronization model must be redone. 

How to define a conventional shower scheme? These are our criteria: i) 
It must be Lorentz covariant and Lorentz invariant, ii) It must be correct 
at LL level. It must be correct at NLL level at least in leading color apprix- 
imation. iii) The soft gluon effect must be correct at least in leading color 
approximation, iv) In every step of the shower the phase space configu- 
ration must be the exact m-body phase space with the exact phase space 
weight. The infrared cutoff parameter should be the only cutoff parameter 
in the algorithm, v) It should "smoothly" work together with the NLO and 
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matrix element matching schemes. 

In this paper we aim to outline some new ideas for a parton shower 
algorithm along our criteria. In the next we will discuss a parton shower 
algorithm based on the Catani-Seymour dipole factorization formulas [3] 
and give a general prescription for the matching of parton shower and fix 
order computation at leading and next-to-leading order level. It was a 
obvious choice to use the dipole factorization because it fulfills our criteria 
and it makes the NLO level matching easier because lots of modern NLO 
fix order computations are based on the dipole subtraction method, for 
example NLOJET++ [5] and MCFM [6]. 

Note, this work is still in progress and we haven't published a proper 
paper about it. Since this is a workshop proceeding there is no way to 
define every details and proof of every statements precisely. 



2. Notation 

2.1. Configuration space 

In order to describe showers, we adopt a vector notation for classical statis- 
tical physics. In order not to confuse classical states with quantum states, 
we use the notation \A) to denote a classical state A instead of using \A), 
which we reserve for a quantum state. The classical inner product is real, 
with (A\B) = (B\A). We denote a state consisting of two incoming partons 
a, b and m outgoing partons with momenta and flavours 

{P, /}a,b,m = {Va.PA, a, ri b p B , b, p U /i, p m , f m } (1) 

by \{p, /}a,b,m)- We can use these as basis states, with normalization 

{{p',f'}*,b,m'\{p, /}a,b,m) = £>mm' <W - rfa.) S b,b> - Vb) 

m (2) 

»=i 

and completeness relation 

1 = E/ /}a,b,m] |{P,/}a,b,m)({P,/}a,b,m| , (3) 

where we introduced the an abbrevation for the integrals and sums, that is 



(4) 
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We can use the basis states to express the differential probability 
M{Pi /}a,b,m) that a general state \A) consists of partons {p, f} a .b,m- 

A({p, /} a ,b,m) = ({p, /}a,b,m|A) (5) 

Then 

l A ) = E / /}^,m] |{p, /}a,b,m)({p, f},,b, m \A) . (6) 

m 

Since ({p, f} m \A) is a probability, a physical state \A) must obey 
the normalization condition (1\A) = 1, where |1) is the vector with 
({f) /}a,b,m|l) = 1- More explicitly the unit vector is 

(!| = E / [ d fa /}».b,m] ({P, /}a,b,m| • (7) 

Given a measurement function F for which the value measured in a 
parton state \{p, /} a ,b,m) is F({p, /} a ,b,m)) it is convenient to define a vector 
l^)by 

(F\{p,f} a ,b,m)=F({p,f} ab!h>m ). (8) 

This is a convenient notation because then the expectation value of F in a 
state A is (F) A = (F\A). 



2.2. Description of color flow 

We will need a description of the color state that is adapted to a description 
of shower evolution. If we use an index notation, to each parton with label 
I there is associated a color index Ai, which takes values 1, . . . , 3 for a quark 
or antiquark and takes values 1, . . . , 8 for a gluon. Thus the matrix element 
has color and spin components and we can write 

Mt:::,t m = E v({c}m) Al - Am M({pj,c} m ) su ..., Sm , (9) 

{c} m 

where the ^({0}™) form a basis for the color space with labels {c} m to be 
explained presently and the M({p, /, c} m ) are expansion coefficients (which 
are still vectors in spin space). In a vector notation, this is 

|M>cb = E l^({ c >™))c ® \M({p, f, c} m )) s . (10) 

{c} m 

Here |A^) c ,s lies in the combined color-spin space while |^({c} m ))c is a 
vector in color space and |.M({p, /, c} m )) s is a vector in spin space. 
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The color basis vectors |V r ({c} m )) c are labelled by a color string config- 
urations {c} m . We define a string configuration to be a set {Si, . . . , S„} of 
one or more strings S. We define strings to be of two types, open strings 
and closed strings. We define an open string to be an ordered set of par- 
ton indices 5* = {h, I2, ■ ■ ■ , l n -i, In}, where l\ is the label of an antiquark, 
l n is the label of a quark, and l 2 , ■ ■ ■ , l n -i are labels of gluons. A closed 
string is an ordered set of parton indices S = {h,h, ■ ■ ■ , ln-i,ln}, where all 
of the indices label gluons and where we treat sets that differ by a cyclic 
permutation of the indices as being the same. 

We describe the color string configuration by assigning a color connec- 
tion label a to every external parton with label i e {a, b, 1, m}. For 
a gluon, Cj = (jr{i), 7r(z)) , for a quark Ci = (0,ir(i)) and for an antiquark 
c, = (7r(i), 0). Here ir(i) is the label of the next parton on the string while 
is the label of the previous parton. This labeling respects that every 
open stings start with a quark and end with an antiquark. The allowed 
string configurations are those in which each parton index is an element of 
one of the strings. 

Now we can define the basis states. We take V^({c} m ) to be a product 

V({c} m ) Al '"- A ™ = V(Si)t A >w V(S 2 ) {A} ™ . . . V(S n ) {A} ™ • (11) 

Here we have denoted the set of color indices represented in string k by 
{A}[ k] = {A h ,. . .,A ln } if string k is S k = {h, ■ ■ ■ ,Z„}. 

We can now define the component factors V(S). We first consider an 
open string. For notational convenience, we suppose that the partons along 
the string are numbered sequentially S = {1,2, ... ,n}. With these labels 
for the partons, we define 

v(S) {A} = [t A n M ■ ■ ■ t^- 1 ]^^ , (12) 

where the t A are the SU(3) generator matrices for the fundamental rep- 
resentation and we take the A\ , A n matrix element of the matrix product 
of the generator matrices. For a closed string with the same parton labels 
(now all gluons) we define 

V{S) {A} = Tr[t Al t A2 ---t A "] . (13) 

A calculation of a cross section involves square of the matrix element. 
Since the basis in the color space is orthogonal up to the leading color terms 
for the matrix element square we have 

|M({p,/k b , m )| 2 = ]T inWa,b, m )| 2 |M({p,/,cK b , m )| 2 + --- , (14) 

{c} a ,b, m 
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where the dots stands for the subleading color contributions. They are 
suppressed by a factor of 1/N%. 

Each term in this sum represents a particular color configuration. Since 
the shower depends on the color configuration, from each term we generates 
a different shower. On the other hand, the leading color approximation is 
useful for nearly collinear parton configurations but is gives a poor descrip- 
tion of the hard part of the event. For better accuracy we should begin 
the showering using the exact matrix element squared reweighted with the 
probability that the hard partons are in a certain color configuration. For 
a given color configuration {c} a ,b,m this probability is 

n f x x_ \v({cK b , m )\ 2 \M({ P j,c} a ^ m )\ 2 

P(|P ' / ' C|a ' b - ) - £ |y( {ckb , m )| 2 |M( {p ,/, c}a , b , m )| 2 • (15) 

{C'}a,b, m 

Now the color configuration assigned matrix element square is 

\M({pJ,c} ath , m )\ 2 = /, C} a>bim ) |.M({p,/}a,b,m)| 2 > ( 16 ) 

where |A4({p, /} a ,b,m)| is the exact matrix element square. 

It is useful to extend the configuration space by adding color informa- 
tion. Introducing the basis vector of the color space |{c} a ,b,m) the new 
configration space is defined as a tensor product that is 

\{P, f, C}a,b,m) = \{P, /}a,b,m) ® |{c}a,b,m) (17) 

and the normalization condition reads 

m 

({c}a,b,m|{c'}a,b,m') = Smm> S Csl c>Jc h c' b Y\.^ c i< ' ( 18 ) 

1=1 

The integral over the state involves a sum over all possible color configura- 
tions 

J [d{p,f,c} a>h , m ]= J [d{p,/} B , b , m ] ■ ( 19 ) 

{c}a,b, m 

The m-parton matrix element square as a vector in the configuration space 
is 

\M m ) - / [d{ P ,f,c} 

a,b,m ] \{P, f, C}a,b,m) \M({p,f,c} a , h , m )\ z , (20) 

where the color assigned amplitude \M({p, f, c} aj b,m)| 2 is defined in 
Eq. (16). 
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2.3. Phase space integrals and Born cross sections 

To define the phase space integral and include the parton distribution func- 
tion it is useful to define an operator in the configuration space. Let us 
define 



T = ^2 [ d &> f> C Kb,m] |{P,/,c} a ,b,m)({P,/,c}a,b,r, 



x /„/A(f|„ Mf) /&/b(%, Mf) r (21) 

x n { (2^p- j+(p?) } (^(w+w-jf-^ft) . 

With this definition we can easy define the phase space integral with the 
proper statistical (1/m!) and flux (l/2rj a r]bPA-PB) factors. Here K is the 
sum of the momenta of the non-QCD particles. 

Now we can define the m-parton Born level cross section as a vector in 
the configuration space 

\a m )=T\M m ) (22) 

where the vector |.M m ) is given in Eq. (20). 

With this notation we can easily define of the Born level cross section 
of any iV-jet quantity. Assuming that F is a iV-jet sensitive quantity then 
its cross section is 

a[F] = (F\a N ) = (F\T\M m ) , (23) 
where (F\ is also a vector in the configuration space 

oo „ 

( F \= S / [ d fe/' C }a,b, m ]F(Ma,b, m )({p,/,c} a , b , m | • (24) 

. AT J 



m=N 



3. Parton shower evolution 

In this section we set up a quite general framework for describing a parton 
shower. We use an evolution variable t that starts at and runs to oo 
and represents someting like log(Q 2 /fc^), where Q is some hard scale. In 
principle, the shower evolves to t — oo but it actually stops at an infrared 
cutoff tf. The evolving shower is represented by a state \A(t)) that begins 
with an initial state |A(to)). The evolution is given by a linear operator 
U(t',t), with 



\A(t)) = U(t,t )\A(t )) 



(25) 
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These operators have the group composition property 

U(t 3 ,t 2 )U(t 2 ,h) = U(t 3 ,t 1 ) . (26) 

The evolution operators preserve probabilities: 

(l\U(t',t)\A) = (l\A) . (27) 

for any state \A). Here (l\A) = 1 if the state is normalized to unit proba- 
bility. 

The class of evolution operators that will use is defined by two operators. 
The first is an infinitesimal generator of evolution or hamiltonian, H(t), 
which can be specified by giving its matrix elements 

({p^C%,b,m'|W(t)|{p,/,c}a,b,m) • (28) 

The second operator is a no-splitting operator N(t', t) that leaves the basis 
states \{p, f, c} ai b,m) unchanged except for multiplying each of them by an 
eigenvalue A({p, /, c} a , b ,m; t' , t): 

N(t',t)\{p,f,c} a , him ) = A({p,/ ) c} B ,b, m ;t / ,t)|{p,/,c} a , b ,m) • (29) 

The evolution operator is expressed in terms of the hamiltonian and the 
no-splitting operators by 

U(t 3 ,h) = N{t 3 ,h)+ [ dt 2 U{t 3 ,h)H{h)N{t 2 ,h) . (30) 

Jt! 

This equation is interpreted as saying that either the system evolves with- 
out splitting from t\ to t 3 , or else it evolves without splitting until an 
intermediate time t 2 , splits at t 2 , and then evolves (possibly with further 
splittings) from t 2 to t 3 . 

There is a relation between the eigenvalues A of N and the matrix 
elements of H. To derive this we apply the evolution operator for a basis 
state and using the normalization condition one can get 

A({p,/,c} a ,b,m;*2,*i) = exp ^-^ dt(l\H{t)\{p,f,c} &ib , m )j • (31) 

Note that the probabilistic meaning of Eq. (31) requires that the exponent 
be negative definite and strictly monotonic in the variable t 2 . 

One of the main purpose of the parton shower programs is to rcsum 
the leading and next-to-leading logarithms which come from the soft and 
collinear regions. In this phase space regions the matrix elements have 
factorization property that any m + 1 parton matrix clement can be written 
as a product of a singular factor and the corresponding m parton matrix 
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element. It is known that the leading and next-to-leading logarithm are 
produced by the 1 — > 2 type splittings. Thus it is enough to consider only 
the 1 — > 2 splittings in the splitting operator H(t) and we have 



This is a very formal definition of the shower algorithm. Let us try to 
define the splitting operators. 

4. Splitting operators 

The formalism outlined above can be used to describe showering with a 
variety of splitting operators. The definition of the splitting operators tht 
we suggest is based on the Catani-Seymour dipole factorization formulas 
and is given by 

({P, f, c}a,b,m+l \H{t) | {p, /, C} a; b,m) 



On the right hand side we have a sum over all the possible emitters (I € 
{a, b, 1, to}). In the general case when we have initial state hadrons the 
incoming partons are included among the allowed emitters. There is a sum 
also over the possible spectator partons (k e {a, b, 1, to} with k ^ I). 
The spitting is parametrized by the splitting variables y, z and <f>. The 
variable y is a virtuality like variable, z is a momentum fraction variable 
and <j) is an azimuthal angle that parametrizes the space-like unit vector k± 
(k\ = —1) which is perpendicular to the emitter and spectator momenta. 

The dimensionless evolution variable t is given by the function 7}^ and 
depends on the {p, /} a ,b,m configuration and on the splitting variables y, 
and z. This function could be any fully infrared sensitive variable. We 
prefer to use the transverse momentum of the emitted partons as defined 
in the centre-of-mass frame of the emitter and spectator. The transverse 
momentum is perpendicular to the momenta of both emitter and spectator. 




(32) 




£ fdz r^6(t + log(T Lk (pi,P k ,z,y)/Q 2 )) 




(33) 
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With this definition, it is boost invariant when we go back to the original 
frame. The strong coupling also calculated at the evolution scale. 

The the splitting kernel is given by the the function S^k{z,y, fi,i)- 
It is derived from the soft and collinear behavior of the tree level matrix 
elements and related to the Altarelli-Parisi splitting functions. In the case 
when we have initial state splitting then according to the backward evolu- 
tion the flavor / ai i = a is the flavor of the new initial state parton. The 
spectator can be either a final state or an initial state parton and of course 
the functional form of the splitting kernel can depend on the type of the 
spectator parton. Since the functional form of the splitting function de- 
pends on the emitter and spectator, one can distinguish four cases. They 
are 



When the emitter or the spectator parton is one of the initial state 
partons, the momentum fraction of the incoming parton changes after the 
splitting. This effect is reflected by the the ratios of the parton distribution 
functions in Eq. (33). The parton distribution functions depend on the 
arbitrary factorization scale y? F . In a shower program \j? F is usually set to 
the evolution variable. However, the exact choice affects only subleading 
logarithms. 

The operator 1Zi t k in Eq. (33) describes the {p, /, c} aj b, TO - * 
{p, /, c} a ,b.m+i transformation and provides the correspondent Jacobians. 

In the dipolc formalism, a splitting always involves two partons. One is 
the emitter (labeled by I in Eq. (33)) that splits and produces two daugther 
partons. The other is the spectator (labeled by k in Eq. (33)) which absorbs 
the recoiled momentum and plays an important role in the color dynamics. 
With this construction, one ensures that the momenta are on-shell and that 
momentum conservation is maintained exactly. Furthermore, the splitting 
probability includes the exact phase space weight. 

How should one choose the spectator parton in a shower picture? In 
keeping with the role of the spectator in the color dynamics for a soft emis- 
sion in the Catani-Seymour scheme, the spectator parton should be color 
connected to the emitter parton. In order to have a good approximation in 
the soft limit we should include the color connection between the emitter 



Si,k(z,y,fi,i,fi,2) = < 



Sfin(z,y,fl,lJl,2) 
Sini(z,y,fl,l,fl,2) 
K Sini(z,y,fl,l,fl,2) 



if Z, k e {1, m} 



if I e {1, ...,m} & k € {a,b} 
if i e {a, b} & fee {1, ...,m} 
if l,k € {a, b} 



(34) 
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and spectator (— TJ • Tk/Tj 2 ). Expanding this operator in l/N c then we 



have 



l _ T 2 = A f, Xi,fc({c}a,b,m) + O 



(35) 



where factor Af l=g = 1/2 for a gluon emitter and for quark or antiquark 
emitter it is Af l=q ^ = 1. The function \i,k is 



Xi,fe({c}a,b,m) — 



1 if Q = (M(Z)) or q = {n'{l),k) 
otherwise 



(36) 



Remember, our color basis are based on open and closed color strings. The 
spectator parton is one of the neighbor parton of the emitter on the color 
string. When the gluon is the emitter there are two possible spectator 
because the gluon always has two neighbors on the color sting while for 
quarks and antiquarks there is always one possible spectator since they are 
always at the end of the color string. 



4.1. Final state splitting with final state spectator 

Let see an example for the splitting functions and for the transformations. 
In the case when both the emitter and the spectator are final state partons 
(I, k e {1, to}) the splitting function is given by 



San(z,y,f 1 ,f 2 ) 



r—u,u,d,d,.. 



6 rfl Srf 2 S q q(z,y) + 5 rfl S g f 2 S qg (z,y) 



(37) 



+ S gfl S rf2 S qg (l - z,y) 

+ 5 9h 5 9h s gg( z ,y) ■ 

Then for the splitting of a quark or antiquark into the same flavor quark 
or antiquark plus a gluon the splitting function is 

2 



S qg (z,y) = C F 



+ 



(38) 



l-z(l-y) 

For the splitting of a gluon into a quark and an antiquark, we define 

S qq (p,z,y)=T R [l-2z(l-z)} . (39) 
Finally for the splitting of a gluon into two gluons the splitting function is 
1 1 



S gg (z,y) = 2C A 



i-s(i-y) i-(i-z)(i-y) 



■2 + z{\-z) 



■ (40) 
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The momentum and flavor transfomation is given by the matrix element 
of the TZ^k, operator and it is 

({P, /}a,b,m+l|^,fe(^,y, /} a ,b,m) 

1 m 



i=i 



8^(p k -(l-y) Pk ) f[ * (4) (Pi- Pi) 



(41) 



i=i 



x 5 (A) (pis - zpi - y(l - z)p k - [2p v p k yz{\ - z)] 1/2 k ± ) 
x <5 (4) (]5;,2 - (1 - z)pi - yzp k + [2p v p k yz{l - z)] 1/2 k ± ) . 

The evolution variable is the magnitude of the transverse momentum, 
that is 

Ti,k(pi,Pk,z,y) = 2pi-p k yz(l - z) . (42) 

4.2. Including heavy quark contributions and SUSY-QCD 

The mass of the top quark and the strongly interacting SUSY particles is 
very big thus the probability of their production is low. In this work we 
don't want to deal with this cases because the phase space requires a more 
complicated treatment. In the future we want to extend our algorithm for 
the heavy and SUSY particles. 



5. Shower cross sections 

In order to calculate the shower cross sections in shower approximation we 
have to define the physical state at the starting (hard) scale to. The hardest 
part of the event is always the simplest process which is kinematically 
possible. For example in e+e~ annihilation the simplest process is e+e~ — > 
qq, or for Higgs production in proton-proton collision it is pp — > H . 

For instance, if we want to study jet production in proton antiproton 
collision then the simplest process is the 2 — > 2 QCD processes and m = 2. 
The to is given by the hard part which could be the the invariant mass of 
the outgoing hard partons. The physical state at the starting scale is the 
simplest state that is kinematically possible 



(43) 
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where |<72) is given in Eq. (22). The shower cross section is just the evolution 
of the initial state 

\a(tf)) =U(t { ,t )\a(t )) , (44) 

where the tf scale is the infrared cutoff scale. And the cross section of any 
jet quantity is given by F is 

<t[F] = (F\D(t f )\a(t f )) = (F\D(t { )U(t u t )\a(t )) , (45) 

where the operator D{tf) represents the hardornization. 

The hadronization is a long distance effect what we cannot calculate 
from QCD but we have some QCD motivated model for it. In principle 
the model is universal (as the hadronization). Once the model has been fit 
to a set of the data {e.g. LEP data) we don't have to retune the model 
for another type of the processes. Although the hadronization model is 
universal, it is strongly coupled with the shower evolution. That means if 
we change the shower evolution we have to retune the hadronization model. 

6. Adjoint splitting operator 

The whole parton shower idea is based on the factorization properties of 
the matrix elements in the soft and collinear regions. In this region an 
m+ 1 parton matrix element can be written as product an m-parton matrix 
element and an universial singular factor. Using this approximtion, one 
can start from a very simple configuration and generate large multiplicity 
multiparton states. 

In our formalism it is essential that the phase space is the exact m-body 
phase space with the exact weight in every step of the shower. Furthermore 
the shower is Lorentz invariant. We saw that applying the operator H(t) 
to an m-parton state, one obtains an to + 1 parton state. The probability 
for the emission is the to + 1 parton phase space weight times a universal 
singular factor. 

Now we can ask if it is possible to do this other way around. Suppose 
that we have generated an to + 1 parton phase space point. Can we define 
an operator that acts on this state and it finds all the possible m-parton 
states with corresponding weights such that these would be the weights of 
an emission if the operator Ti was applied on these m-parton states. We 
claim it is possible to define this operator since the factorization works 
exactly this way. 
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Let us define the adjoint splitting operator T0{t) that obeys to the 
following condition: 

(F\H(t)T\A) = (A\H\t)T\F) , (46) 

where the states \A) and \F) are arbitrary and the operator T makes the 
integrals phase space integrals as it is defined in Eq. (21). We can see im- 
mediately that the operator TO (t) always decreases the number of partons 
because TL(t) always increases it. Since in Tt(t) we consider only 1 — > 2 
splittings in TO {t) we have only 2 — > 1 merging. Thus there is only one 
non-zero matrix element 
(l|Kt(t)|{p, /, c} a , b , m ) 

f - (47) 

= / [d{P,/,c} a ,b,m-l] ({p,/,c}a,b,m-l|W t (t)|{p,/,c} a ,b,m) • 

From Eq. (46) we can immediately see that for several emission we can 
write 

(F\H(t m )H(t m - 1 )---H(t 3 )r\A) = (A\TL\h)---TL\t m _ 1 )TL\t m )T\F) . 

(48) 

Using Eq. (33) and from the Catani-Seymour dipole subtraction terms 
[3] for the non-zero matrix elements of operator 7~0(t) we have 

({P, f, c}a,b,m|H t (t)|{p, /, c} ajb ,m+l) 
1,3 k&rf J 

pairs ^ ' 

x A fij xtj,k({c}^b,m) S(t + \og(T ll . k (p l ,p :j ,p k )/Q 2 )) 

X ({p, /, c} a , b ,m | Qij.k | {P, /, c} a ,b ,m+l J 

Here we have sum first over all the possible pairs {i, j} of partons that 
could be the daughter partons from the splitting. The daughter partons 
must be color connected. Then we sum over possible choices k of the 
spectator parton. The Dirac delta function maps out the dimensionlcss 
evolution variable t, which is the logarithm of the transverse momentum 
Tij,k (Pi j Pj 7 Pk ) in the splitting. We have the spin averaged splitting function 
(Vij t k) that depends on the momenta of partons i,j, k G {a, b, 1, ..m + 1} 
and the flavors of partons i and j. It contains a factor of the strong coupling 
a s (kj_). The color connection between the emitter and spectator (— Tj 



T)-/Tfj) is considered only at leading color level, that is Af i:j Xij,k- The 
momentum, flavor and color transformation prescription is specified by the 
operator Qi^k- 



February 2, 2008 5:1 Proceedings Trim Size: 9in x 6in 



Ringberg05 



15 

7. Matching parton shower and Born level matrix elements 

In Section 5 we defined the shower cross section as the shower evolution of 
the kincmatically simplest hard process, for example the 2^2 processes 
in hadron-hadron collisions. It is obvious that this approximation is very 
poor if we are interested in the cross sections other than the 2-jet cross 
section. In order to have a better approximation one should consider higher 
multiplicity exact matrix elements. 

Let us start with the shower cross section expand the first n steps of 
the evolution operator 

|<7(t f )) =N(t { ,t 2 )\a 2 ) 

+ Y, / dt * dt 4--- dt m N{t u t m )H{t m ) 

m = 3 Jt 2 Jt 3 Jt m -! 

x N(t m , t m -i)W(t m -i) • • • N(U, t 3 )H{t 3 ) N(t 3 ,t 2 )\a 2 ) (50) 
+ f 'dts /" dt 4 ■ ■ ■ f 1 dt n U(t { , t n )H(t n ) 

x N(t n , t n _i)W(t„_i) • • • N(U, t 3 )H{t 3 ) N(t 3 ,t 2 )\a 2 ) . 

In every term we have a string of operator products H(t)N(t,t'), which 
provide the probability of non-emission between "times" t' and t followed 
by a splitting at t. This corresponds to the approximate matrix element 
squared times the proper Sudakov factors. We want to replace the approx- 
imate matrix element with the exact m-parton tree level matrix element 
while keeping the proper Sudakov weights. 

We use the approximation of the matrix elements in the soft and 
collinear regions when emissions are ordered. One can show if the state 
\{Pi /) c }a,b,m) was generated according to the shower procedure then the 
following is a good approximation 

(■Mm|{P,/,c} a ,b,m) ~ (An(*f , t 2 ) \{p, f, c} ai b, m ) 

ee [ u dt 3 [ u du ■■■ T dt m f[ a ;Jf R \. (51) 

x (M 2 \H\t 3 )HHt 4 )---HHt m )\{p,f,c} a ^ m ) . 

It is clear that knowing the state \{p, f, c} ai b, m ), the approximate matrix 
element (A m (tf,t 2 )\{p,f,c} a ^ iTn ) is calculable according to the definition 
of operator TV (t). However this could be a very big sum with lots of terms. 
One can immediately see that this approximation is valid only in those 
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regions where the emissions are ordered and the phase space is dominated 
by soft and collinear radiation. 

In the next step we reweight every term in Eq. (50) by a matrix element 
factor. Whenever (A^ m | is known and the state \{p, /, c} a ^,m) is generated 
by the shower procedure, the following ratio is always well defined: 

i r f \ ( i \ {M m \{p,f,c} a , h ,m) , . 

w M ({p,f,c} a , h , m ,ti,t 2 ) = ' , r- . (52) 

(Am(t{,t2)\{p, J,c) a: b,mj 

When \M m ) is unknown, we can simply set wm{{p, /, c} a ,b,m, *f , t 2 ) = 1. 
Let us define the matrix element reweighting operator 

V^M^f,^) = / /, c}a,b.m] \{p, f, c}a,b.m)({p, /, c}a,b,m| 

m J (53) 

X tfM({P,/,c} a)b)Tn ,tf,t 2 ) • 

Acting on an m-parton state generated by the shower algorithm, this oper- 
ator replaces the approximate squared matrix element by the exact one (if 
m is small enough that the squared matrix element is known). Acting on 
a state with larger m it acts as the unit operator. Using this reweighting 
operator, the improved shower cross section is 

\<J M (t{)) = N{t{,t 2 )\a 2 ) 

+ H/ d<3 / dt i--- dt m N(t f ,t m )W M (tf,t 2 ) 

771 = 3 * 2 * 3 Jtm-1 

x WCOMWm-i) ' ■■H(t 3 )N(t 3 ,t 2 )\a 2 ) (54) 
+ [ dta [ dU---( dt n U(t { ,t n )W M (t { ,t 2 ) 

Jt 2 Jt 3 Jt n -! 

x H(t n )N(t n ,t n ^) ■ ■ - H(h) N(t 3 ,t 2 )\<i 2 ) . 

Here it is assumed that the matrix elements lA^), • • • , \M n ) are known. 
The nice feature of this formulae that the event is generated according to a 
shower algorithm starting form the trivial 2^2 kinematics and whenever 
the exact m-parton matrix element is known the event is rewieghted by the 
ratio of the this tree level amplitude and the approximate shower amplitude. 
This form is suitable for numerical implementation. 

With the exact matrix element correction we don't change the leading 
and next-to-leading logarithms of the cross section. It is easy to see if we 
rewrite Eq. (54) in an equivalent form by adding and subtracting the same 
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terms. Thus we have 



<Jm{M)) = U(t { ,t 2 )\a 2 ) 

+ Y1 / d< 3 / dt A --- dt m U(t f ,t m )[W M (tf,t 2 ),H(t m )] (55) 
x N(t m , t m -i)H(t m -i) ■ ■■H(t 3 )N(t 3 ,t 2 )\a 2 ) 



[W M {tf,t 2 ), H(t m )] = W M (tf,t 2 )H{t m ) - H(t m )W M (tut 2 ) ■ (56) 



In Eq. (55) in the terms with higher multiplicity matrix element correction 
the leading (LL) and next-to-leading logarithms (NLL) are removed by the 
commutator leaving only subleading logarithms and finite pieces. All the 
LL and NLL logarithms are accumulated in the the first term which is 
actually the standard shower. 

7.1. Matching with Sudakov reweighting 

In this subsection we rewrite our matching formulae in an equivalent form. 
This form is less efficient from the point of the numerical implementation 
but it is very useful for further studies. 

After some algebraic manipulation one can find that Eq. (54) can be 
written in the following form: 

| CM (if)) = \(TA(tf)) 



Jt 2 

where the state \a m ) is defined in Eq. (22) and the Sudakov reweighting 
operator is 



and the commutator is defined in the usual way, that is 




N(t u t 2 )\a 2 ) + / 

dt m N(t{, t m )W&{tu t m , t 2 ) \<T m ) 



(57) 




m 




X 



(M 2 \N(t; 



■3, t 2 )HHt 3 ) ■ ■ ■ N(t, f m -i)ttt(f)|{ p) /, c} a ,b,m) 
(Am(tt, t 2 )\{p, f, c} aj b,m) 



(58) 
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This operator provides all the possible emission history with the proper 
Sudakov factor. 

7.2. Slicing method 

The slicing method was originally defined by Catani, Krauss, Kuhn and 
Webber [4] for e + e~ annihilation and was later was generalized for pro- 
cesses with initial state partons by Krauss [7] and Mrenna and Richardson 
[8] . Nowadays this method is very popular and one can find several imple- 
mentations in the literature [9]. It is usually called the CKKW-method. 
The basic idea is to divide the evolution region two parts, a hard part and 
a soft part. In the hard region one uses the exact matrix element sup- 
plemented by an appropriate Sudakov factor. The soft part is dominated 
by soft and collinear radiation so the shower cross section provides a good 
approximation. 

Defining the matching scale if > i; n i > io and using the group decom- 
position property of the evolution operator then we have 



Now, we can apply the matching formulae for t/(ij n i,i2)| 02(^2)) according 
to Eq. (57). The CKKW-method uses a further approximation for the Su- 
dakov reweighting operator. Instead of using (57), the emission history and 
the Sudakov factors are defined according to the k± jet finding algorithm. 
With this method, a unique emission history is associated with each given 
n-parton configuration. For the CKKW cross section we have 

|fCKKw(if)) = U{t i ,t ini )N(ti n i,t 2 )\<J2) 



The CKKW method has the advantage that the Sudakov reweighting 
is very simple. Thus the method is easy to implement. However, we have 
to pay the price that the result depends on the arbitrary matching scale. 
Although the t; n i dependence is canceled at next-to-leading logarithmic 
level, this dependence could still could be a source of uncertainty [8] . 



a(t { )) = U(t f ,t ini )U(t ini ,t 2 )\a 2 (t2)) 



(59) 




n—l r ti 



(60) 
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8. Matching parton shower and matrix elements at 
next-to-leading order level 

Matching the parton shower and the NLO computation is a long outstand- 
ing problem. Recently some very important progress has been made. One 
example is the program of Frixione, Nason, and Webber [10], which so far 
has been applied to cases with massless incoming partons but not to cases 
with massless final state partons. The other example is that of [11], which 
concerns three-jet observables in electron-positron annihilation. 

In the next we review a general algorithm for matching fix order NLO 
computation and parton shower. 

8.1. NLO cross section 

Before we discuss the matching at NLO level it is useful to have a brief 
review of the structure of the NLO fix order calculation. According to the 
dipolc method in the general case the cross section is 



The first term da B is the Born contribution. The second term is an N + 1 
parton phase space integral which gives the contribution of the real radi- 
ations da R . Since da R is singular in some region of the phase space and 
this singularities are regularized by the dipole subtraction terms da A . The 
dipole contributions are based on the dipole factorization formulas [3]. 

The next term is the contribution of the virtual graphs da v with N 
outgoing external partons. They have also infrared divergences in d = 4 
dimension in terms of 1/e (where epsilon is the parameter of the dimensional 
regularization) but these divergences are cancelled by da B Furthermore 
we have some finite ./V-parton contributions da B <g> [K + P(fip)]. 

8.2. Cross section with shower 

Let us calculate the cross section of an infrared safe -/V-jct quantity. Assum- 
ing we know \<tn) and | a jv+i) the matrix element improved cross section 




(61) 
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from Eq. (57) is 




(F N \D(t { )\a A {t{)) = / dt N (F N \D(tf)N(tf,t N )W A (ti,t N ,t 2 )\a N ) 



where operator D{t{) represents the hadronization and function Fpf defines 
the iV-jet observable. Since Fn is an infrared safe quantity it is safe to ex- 
pand this formulae in a s and drop the power corrections which are provided 
by the hadronization. After the expansion one can show that 



(F N \D(t f )\a A ) = [ da B (l + C+^W A 1) )+ [ [da R - da A ] 

JN V 27T / Jn+1 

+ 0(a 2 s )+0(lGcV/^) ■ 



On the right hand side we have two integral. The second term is the 
N + 1 parton integral of the real radiation like in Eq. (61). The first term 
consists of three contributions. We got the Born term da B . There are some 
NLO contribution da B W A 1] that comes from the expansion of the Sudakov 
reweighting operator W A . 

The contribution da B C is the error term of the leading color approxi- 
mation. Since the splitting operator Ti(t) is valid only at leading color level 
in the soft gluon limit this error term generally is not zero. It is a leading 
order contribution but it is only sublcading color contribution 



On the other hand only next-to-leading logarithms contributes to this term. 
One can neglect this term but with minor redefinition of Eq. (62) it is 
possible to manage that this term vanishes [12]. In this proceeding we 
simply neglect this term. 

One can immediately see that the expression on the right hand side of 
Eq. (63) is almost the NLO cross section. We got the Born term and the real 
contributions right but we missed the virtual and some finite NLO contri- 
butions. In order to match the shower and the NLO fix order computation 



(62) 




(64) 
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we have to modify Eq. (62). Thus we have 

(F N \D(tf)\trNL0(tf)) = f dt N (F N \D(t i )N(t i ,t N )W A (tutN,t2)\(JN) 
Jt 2 

,t t 

+ / dt N+1 (F N \D(ti)U(ti,t N+1 )W A (t{,t N+1 ,t2)\(TN+i) 
Jt 2 

ftt 

+ / dt N (F N \D(t { )U(t { ,t N )W A {t u t Nl t 2 )\a% LO ) . 

(65) 

Here we have an extra term that is basically the shower evolution of the 
state \<Jn L °)- This state provides all the missing terms for the NLO match- 
ing 

[4c}a,b,iv]« Z ' |{p,/,c} a , b , J v) = \M N \ 2 

2n (66) 
+ \M N \ 2 ®(K + Pfa)) + \\M N \\_ loov + \M N \ 2 ® 1(e)] 

where tv | 2 and \M.n\\-i oop are the tree and 1-loop level matrix element 
squares, respectively and the insertions operators /(e), K and P(hf) are 
defined according to Eq. (61). 

However this matching formulae is very simple and general we still have 
some caveats. Even if the problem with the error term of leading color 
approximation is solved the matching is exact only for the 2^2 like 
processes. That is because in the shower the emissions are ordered and we 
use only the leading order splitting kernels. This ordering cut out a certain 
region of the available phase space which is usually considered in the fix 
order NLO calculations. On the other hand this error is negligible because 
in this region the NLO computation is not reliable. 



8.3. Slicing method at NLO level 

In order to have an simpler Sudakov reweighting procedure one can do 
the NLO matching with slicing method. This algorithm is completely 
worked out for processes in e+e~ annihilation and the details can be find 
in Rcf. [12]. 

9. Conclusions 

In this paper we outlined a new parton shower algorithm based on the 
Catani-Seymour dipole factorization formulas. We also introduced a new 
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formalism for describing the shower and we find it very powerful. The 
main advantages of the this new algorithm arc the followings: i) It is 
Lorentz invariant and covariant formalism, ii) The evolution is a trans- 
verse momentum ordered algorithm and it is managed that the partons are 
always onshell, the momentum conservation is maintained and the phase 
space weight is the exact weight in steps of the evolution, iii) There are 
no ambiguous technical parameters. The only parameter is the infrared 
cutoff scale parameter, iv) Improved soft gluon treatment. Although the 
soft gluon limit is still considered in large N c limit we could make some 
improvements by taking into account not only the final-final but also the 
initial-final and initial-initial state color connections. 

We derived a more general formulae for matching Born level matrix 
elements and parton shower. Our matching algorithm is free form any 
arbitrary matching parameter but we show that the CKKW method is a 
specialization of our matching procedure. 

Since the shower evolution is based on the dipole factorization, the 
matching of the shower and NLO fix order computation (based on the 
Catani-Seymour subtraction method) is rather straightforward. 
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